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We examine the non-Markovian dynamics in a multipartite system of two initially correlated atomic qubits, 
each located in a single-mode leaky cavity and interacting with its own bosonic reservoir. We show the dom- 
inance of non-Markovian features, as quantified by the difference in fidelity of the evolved system with its 
density matrix at an earlier time, in three specific two-qubit partitions associated with the cavity-cavity and 
atom-reservoir density matrices within the same subsystem, and the cavity-reservoir reduced matrix across the 
two subsystems. The non-Markovianity in the cavity-cavity subsystem is seen to be optimized in the vicinity 
of the exceptional point. The CHSH-Bell inequality computed for various two-qubit partitions show that high 
non-locality present in a specific subsystem appears in conjunction with enhanced non-Markovian dynamics in 
adjacent subsystems. This is in contrast to the matching existence of non-locality and quantum correlations in 
regions spanned by time t and the cavity decay rate, A c for select partitions. We discuss the applicability of these 
results to photosynthetic systems. 

PACS numbers: 03.65.Yz, 03.65.Ta, 03.67.Mn, 42.50.Lc, 71.35.-y 



I. INTRODUCTION 

One of the most intriguing feature that underpins the many promising applications |Q]] of quantum theory is that it cannot be 
reduced to local realistic models in which the outcomes of local measurements are determined in advance by hidden variables 
J2]. Non-locality which is linked to the violation of Bell inequality of any form, is considered to arise from complementarity or 
the impossibility of simultaneous joint measurements of observables. The Clauser-Horne-Shimony-Holt (CHSH) inequality 
is the only extremal Bell inequality with two settings and two outcomes per site, and involves the comparison of predictions of 
quantum theories with theories rooted in local realism. 

Entanglement and non-locality are distinct entities which are based on different frameworks. While entanglement is invariably 
linked to the mathematical formalism of the abstract Hilbert space, non-locality is founded on correlations of variables that, while 
amenable to experimental verification, has rather debatable origins due to the paradoxical, non-Boolean nature of quantum logic, 
with striking inconsistencies with the deterministic structures that underpins classical logic. The non-objectivity-non-locality 
issue invariably leads to inherent difficulties in formulating a rigorous definition of non-locality. Accordingly, the violation of 
any Bell inequality may imply the presence of non-locality, however a non violation does not guarantee the local or nonlocal 
nature of the correlations under study. Interestingly, the violation of a Bell inequality does not imply the existence of non-locality 
within the system where Bell test was performed, but it is a requirement |4] that non-locality is observed in the vicinity of the 
quantum system under observation. While non-locality is a feature which reflects entanglement, the converse is not true i.e., 
not all entangled states reveal non-local features. In other words, non-locality is a sufficient feature but not necessary to reveal 
entanglement ^ . 

A closely associated problem that has yet to be examined in depth is the role of non-divisible quantum maps, first proposed 
by Sudarshan et al, Jal and examined in other related works 01 17 UT2I1 . on the violation of the CHSH inequality and the 
association between non-locality and non-Markovianity. The challenges of relating non-Markovian dynamics and maps that 
are not completely positive (CP) and which describe the evolution of open quantum systems is well known fl3ll . The quantum 
Markov process as opposed to its non-Markovian counterpart is described by a subset of completely positive intermediate time 
dynamical maps ifTitl . is limited by its strict applicability to weak system-environment interactions. In a dissipative environment, 
as is the case in this study, non-Hermitian effects are introduced and possibly a non-Markovian evolution involving some back 
flow of information to the system ensues. It would be interesting to examine non-Markovian processes of open quantum systems 
with non-divisible maps and which incorporate memory feedback mechanisms within the realm of non-locality. In a recent work 
lfl5tl . the link between non-Markovian revivals and non-locality for a qubit-oscillator system was examined with results showing 



'Electronic address: thilaphys@gmail.com 
^Electronic address: ushathirthahalli@gmail.com 



2 



a strong correlation between revivals of non-locality and non-Markovian evolution dynamics of the studied system, nevertheless 
these trends were noted as not being indicative of a non-Markovian basis to changes during violation of the CHSH inequality. In 
an earlier work fl^l . the association between entanglement and non-locality was examined from a dynamical perspective, and 
closed relations were obtained for independent environments. 

The main interest in non-Markovian processes arise from the fact that Markovian dynamics is typically only an approxima- 
tion, which is no longer valid when considering shorter time-scales and/or stronger system-environment couplings as in light- 
harvesting systems where non-Markovianity appears significant under physiological conditions lfl7i [Till . The non-Markovian 
measure has several interpretations with its evaluation based on a variety of physical measures depending on the evolution dy- 
namics of the system under study. One popular computable measure was proposed by Breuer and coworkers fl9ll . where Marko- 
vianity is employed as a characteristic of the dynamical map p(0) h> p(t) = <t>(f, 0)p(0), and is associated with the decrease in 
the trace-distance which quantifies the distinguishability, D[pi,pi\ between two system states, p\, pi ■ The distinguishability 
measure does not increase under all completely positive, trace preserving maps, hence <x = dD[p\,p2, t]/dt is negative (positive) 
when information flows from (to) the system to (from) its environment. Consequently, the increase of trace distance during any 
time intervals is taken as a signature of the emergence of non-Markovianity over the period of dynamical evolution, and has been 
employed to examine dynamics of interactions in light-harvesting systems lfl8ll . The quantitative estimate of non-Markovianity 
is based on the cumulant of the positive information flux involving a maximization of all pairs of initial states to determine the 
largest amount of information that can be recovered from the environment lfl9tl . Wolf et al. |f20[l introduced an alternative mea- 
sure of non-Markovianity based on divisibility, and subsequently Rivas et al. UM also proposed a measure of non-Markovianity 
based on deviations from divisibility and the characteristics of quantum correlations of the ancilla component of an entangled 
system evolving under a trace preserving completely positive quantum channel. The subject of equivalence of the two measures 
of non-Markovianity for open two-level systems has been topic of much interest in recent works I22II23I1 . with the results indi- 
cating a dependence on the generic features of the quantum system for the two measures to be reconciled. It is to be noted that 
both measures are based on deviations from the continuous, memoryless, completely positive semi-group feature of Markovian 
evolution as is the case in other measures introduced recently. 

Lu et al. l24ll defined another measure of non-Markovianity using the quantum fisher information flow which quantifies the 
computable phase of the qubit undergoing decoherence due to its environment. Recently the fidelity of teleportation l25ll was 
seen to exhibit increased oscillation of fidelities in the non-Markovian regime of photosynthetic systems. Revivals of quantum 
correlations can exhibit revivals even in absence of back-action from the environment to the system, this feature was linked to 
the non-Markovian character of the map using a suitable quantifier of non-Markovianity in Ref. l26ll . Rajagopal and coworkers 
rf28tl have explored the generality of the Kraus representation for both the Markovian and non-Markovian quantum evolution, 
highlighting the role of the fidelity measure F[p(t),p(t + r)] l29ll as being able to track salient properties of an evolved density 
matrix p{t+f) as compared to a reference density matrix p(t). This direct approach offers convenience in gauging non-Markovian 
signatures based on the fidelity difference F[p{t),p(t + t)] - F[p(Q),p(T)]. 

Our aim is twofold: (a) to examine the general relationships between non-Markovianity, non-locality and non-classical corre- 
lations l27ll in a multipartite system consisting of a qubit interacting with a structured reservoir in the form of a cavity in contact 
with its own reservoir. This multipartite arrangement provides different channels of decoherence and dissipation: one that occurs 
between the qubit and cavity and the second, between the cavity and reservoir, (b) to investigate the characteristic behavior of 
the entities in the vicinity of the exceptional points which are topological defects ll30ll . which occur when two eigenvalues of 
an operator coalesce when selected system parameters are altered. As a result, two mutually orthogonal states merge into one 
self-orthogonal state, resulting in a singularity in the spectrum [f30h - Using the qubit-cavity -reservoir multipartite arrangement, 
we seek to examine the links between non-Markovianity, violation of the CHSH-Bell inequality and non-classical correlations 
in regions governed by the system parameters within various two-qubit partitions. 

This paper is organized as follows. In Section [II] we derive analytical expressions for the atom-cavity -reservoir model with 
cavity losses based on a phenomenological master equation and the quantum trajectory approach. In Section lHAl the entangle- 
ment dynamics of a multipartite system of two noninteracting subsystems, each consisting of the two-level atom coupled to the 
cavity which interacts with its reservoir source is examined, with expressions obtained for several bipartite two-qubit density ma- 
trices. In SectionUm analytical expressions of fidelities for several two-qubit partitions near the exceptional point are provided. 
Analysis of qualitative features of non-Markovianity based on the fidelity difference measure of specific subsystem are analysed 
in Section [TV] In Subsection II V Al qualitative features of non-Markovianity of the cavity-cavity subsystem obtained in Section 
HVlis compared with those obtained using the trace-distance difference measure. In Section [V] we make observations based 
on the non-Markovian results in Section HVl and results of non-locality computed using the CHSH-Bell inequality. In Section 
IV Al we compare results of Bell non-locality and non-classical correlations for given partitions, and discuss the applicability of 
the results obtained in this study to photosynthetic systems in Section[VTJ Finally, Section [VHl provides our summary and some 
conclusions. 
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II. ATOM-CAVITY-RESERVOIR MODEL WITH CAVITY LOSSES 

We consider a multipartite system of Hamiltonian Hj which consists of a two-level atom coupled to a single-mode leaky 
cavity s, which in turn is interacting with its own source of bosonic reservoir r {ti=\) 

H T = H s + H r + H, (1) 

H s = cl>q cr + cr_ + u> c a ' a + V(cr_ a ' + cr + a) (2) 

N 

H r = 2^ Ukb\bk, Hi = ^(<fkab" k + <P* k bka T ), (3) 
k k=\ 

where ojq is the atomic resonance frequency, <x + (cr_) denotes the raising (lowering) operator of the atom, and a (a 1 ) annihilates 
(creates) a photon with frequency u> c in the cavity mode. The operator bk (bf) annihilates (creates) a photon with frequency 
ojk in k-ih mode of the reservoir. V is the coupling constant between qubit and cavity and <pk is the linear coupling between the 
cavity and reservoir. 

The density operator p of the quantum system associated with the total Hamiltonian, Hj (Eq. (Q]i) is obtained using the 
generalized Liouville-von Neumann equation ^ = —i£p, where X generates the map from the initial to final density operators 
via a Liouville superoperator <D(f, 0): p(0) i-> p(f) = <D(f, 0)p(0). For the simple reversible and non-dissipative situation, and one 
in which the state commences in an uncorrected combination of the atom-cavity subsystem and reservoir, the reduced dynamics 
of the atom-cavity subsystem can be obtained in the Kraus representation, by taking the trace over the reservoir degrees of 
freedom 

p(/)=^^(/)p(0)^(0, (4) 

i 

where the operators K, constitute the Hilbert space of the subsystem, with 2/ K] K, = I, and / is the identity matrix. The map 
in Eq. (0]i projects a density operator to another positive density operator, and hence is a positive map. The separation of the 
total quantum system into a subsystem of immediate interest, and the reservoir environment which is generally considered to 
be in thermal equilibrium forms the basis of Eq. This is not the case in most practical situations when both the Liouville 
operator and the map <t>(f, 0) (for negative f) are not well defined as the dynamical evolution of the quantum state of the total 
system (atom, cavity and reservoir degrees of freedom) becomes irreversible and dissipative in nature, making it cumbersome 
for realistic modelling of even the atom-cavity subsystem. 

Using the projection-operator partitioning method proposed by Feshbach l3lll . the total Hilbert space of Hj (Eq. ((TJ) is divided 
into two orthogonal subspaces generated by the projection operator, V and its complementary projection operator Q — 1 - P. 
This allows the convenient study of the subspace of interest (atom-cavity in our case) within the total quantum system specified 
in Eq. (OJ. For our system, V — cr + <X- + a'a and Q = 2a bJ>k, such that fQ = QP=0. This means that the operators of 
the atom-cavity subsystem and the reservoir subsystem (which acts as the dissipative source) commute. The density operator 
associated with the qubit-cavity system is obtained using 

p s (t) = PpP, (5) 

where p on the RHS of Eq. © is the density operator of the total system. An excitation present in the two-level atom a (cavity 
c) is denoted by |l) a (|l) c ), whilst |0) a (|0) c ) correspond to the ground state atom (vacuum state of the cavity). The subspace V 
is spanned by (|l) a , |0) fl , |l) c , |0) c ) while Q is spanned by (|0),. , |l) r ). The reservoir state |0) r =I~[i |0)* denotes the vacuum of the 
reservoir, however its ortho-complement state, |l) r needs further explanation. We note that other than evolving under the action 
of the Hamiltonian associated with the atom-cavity subsystem (H s , Eq. (fJI), the reservoir states also evolve under the reservoir 
Hamiltonian, H r in Eq. OJ. Hence we define the collective state of the reservoir sink, |l) r , as a series of superposed terms, 
consisting of the state with a single excitation, |l) r .= 2* <Pk\lk)r where Rq = VSt \<Pk\ 2 , and other states orthogonal to it, |l) r . 
(i = 2,3...) as follows 

n>,=- y>< n> ri (6) 

where \i is the probability amplitude associated with the existence of an excitation in the collective reservoir state, with the atom 
and cavity, both remaining in the zero state. |l*) r denotes an excited oscillator in the A:-th reservoir mode, with oscillators in all 
other modes remaining in the unexcited state. 

Instead of a full microscopic derivation of the system dynamics associated with the Hamiltonian, H T (Eq. ([TJ), we incorporate 
cavity losses by means of a phenomenological master equation which provides reliable results when the spectrum of the reservoir 
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is approximately flat, as justified in a recent work (32j which compared this model and a microscopic system-reservoir interaction 
model. The phenomenological master equation involving the reduced density matrix, p, appear as 

^-p s = i(p s H s - H s p s ) + ^■(2cT-p s o-+ - lt+ct-Ps - p. s cr +C r_) + ^■(2ap s a i - a'ap s - p s a'a), (7) 
at 2 2 

where A a (A c )represents the atomic (cavity) decay rate. The master equation in Eq. © can be further simplified (setting A a =0) 
as follows 

^p s = i(p s H' s - H' s p s ) + A c ap s a\ (8) 

H[ = H s - fl'fl (9) 

where the second term in the Hamiltonian H' s constitutes a decay component which contributes to the non-Hermitian features 
in the system. The last term in Eq. (O denotes the action of the jump operator which forms the basis of the quantum trajectory 
approach l33l - l36ll . In this approach, incoherent processes associated with non-Hermitian terms are incorporated as random 
quantum jumps which results in the collapse of the wavefunction. The net effect is the transfer of states from one subspace to 
the other, P — > Q. The density operator is obtained by taking an ensemble average of a range of conditioned operators l33tl 
at the select time t. It is important to note that while Eq. ^ describing the dynamics of the atom-cavity system is Markovian 
within one subspace, non-Markovian dynamics may arise in the interaction between the atom (cavity) and reservoir in different 
subspaces. 

In order to obtain analytical solutions, we rewrite H\ in Eq. (0 (c.f. H s , Eq. (f2]i) as 

A C .L .[. 

H s = too c+ <T- + (to c + 6- i — )a'a + V(cr_ a' + <x + a) , (10) 
<*c f \<Pk\ 2 

8- i—= g{u k )dco k (11) 

2 J to — to k 

where 5 arises due to renormalization and energy integration is based on the assumption that the reservoir state energies are 
closely spaced. We intentionally ignore detailed parameters, such as spectral density and temperature, of the reservoir system 
for simplicity in analysis of the model under study. These boson bath parameters invariably determine the decay rate A c as is 
obvious in Eq. (0. Moreover information related to the distribution of the phonon bath frequency including the cutoff frequency 
will influence the bath memory time, and accordingly the decay rate A c incorporates leakages between subspaces. The latter 
processes are expected to result in a complicated mix of the distinct subspaces of f and Q, with bearings on non-trivial non- 
Markovian dynamics between two subsystems present in the Hilbert space of Hj (Eq. ([T]i), this will be revealed in forthcoming 
Sections. 

Considering the resonant condition coo - co c + 5 and the presence of a single initial excitation in each subsystem, we write the 
multipartite state of the atom-cavity-reservoir as 

lfc> = U> fl ®|0>c®|0> r , (12) 



Following the quantum trajectory dynamics model 113311 . Eq. (fT2l evolves under the action of the Hamiltonian in Eq. (HJ as 
follows 

tyt) = 6U>« I0) c |0>, + rit |0>J1> C |0> r + X t |0>, |0) c |l) r , (13) 

where \^,\ 2 {\tj,\ 2 ) is the probability that the excitation is present in the atom (cavity). The probability amplitude, Xu introduced 
earlier in Eq. © is evaluated using \x,\ 2 = 1 - |£ f | 2 - \t],\ 2 - 

For the initial condition, £o=l> ?7o=0, the analytical expressions for and i], can be obtained using 

^§ = -iVm (14) 
at 

dm .... A c 

-dT = - ,Vfi " 2 th 

To simplify further analysis of the problem, we set |£,| 2 = 1 — p, \ri,\ 2 = q and \x t \ 2 = p - q. Here q(t) is the probability of the 
atomic qubit exchanging quantum information with the cavity at time f, and 1 - p(t) is the probability that the atomic qubit will 
remain in its excited state, provided it was in the excited state at t =0. The exchange mechanism obviously includes a dissipative 
measure, jd=p-q which is non-zero in the event that the reservoir states coupled to the cavity gets excited. Accordingly, there 
exists a probability p — q that the atomic qubit decays without exciting the cavity states, and we note that in the presence of 
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non-Hermitian exchanges p > q, yd >0. In general p(t) grows with time, and as will be shown later, it can incorporate feedback 
mechanisms associated with the non-Markovian dynamics of the multipartite state. 
Eq. ([Pil l yields the analytical forms for p, q and appear as 



p - 1 - e 



-Act/2 



cos Qf H sin Qf 

4Q 



q = e -<W2Zl sin 2 Qf5 



(15) 
(16) 



where the Rabi frequency 2Q = (4V 2 - (— ) 2 )'^ 2 - Similar expressions as in Eq. (15\ c an be obtained via inversion of a Green's 
function as detailed in an earlier work on the dissipative two-level dimer model 0711 . We note the existence of two regimes, 
depending on the relation between V and Ac- 1 
regime, and at the exceptional point we obtain 



depending on the relation between V and A c . The range where V > 4f (V < 4r) applies to the coherent (incoherent) tunneling 



Q 



A c 

0, V = — 
4 



At this point , both the coherent and incoherent tunneling regime merge and we obtain 



(17) 

(18) 
(19) 



The exceptional point is a topological defect which is present in the vicinity of a level repulsion BOH and unlike degenerate 
points, only one eigenfunction exists at the exceptional point due to the merging of two eigenvalues. The critical temperatures at 
which exceptional points occurs in photosynthetic systems was identified recently 12511 , In the restricted subspace subspaces of 
V occupied by the atom-cavity system, we note that the exceptional point in Eq. ([T7| i may take on a range of values as the cavity 
decay rate, A c is linked to almost continuous range of reservoir attributes. 



A. Entanglement dynamics of a multipartite system incorporating non-Hermitian terms 



Here we extend the composite system in Eqs. (fT3l l to examine the entanglement dynamics in a system of two noninteracting 
subsystems, each consisting of the two-level atom coupled to the cavity which interacts with its reservoir source. To simplify 
the analysis, we consider that the subsystems have identical environment with respect to their reservoir characteristics and cavity 
decay attributes, and assume a multipartite system in the initial state 



m = (a \\ ai % 2 ) + b |0 Ol U>) |0 C1 C2 > |0 n 0„) , 



(20) 



where a, b are complex parameters, and the cavity c and reservoir r are present in their vacuum state. The two subsystems are 
labeled as 1 or 2. 

The composite system, evolves as 



a [f ( |l> fl| |0> C] |0> n + m |0> flI |1> C1 |0),, + Xl |0> ai |0> C1 |l> n ] |0 a2 C2 rz > 
+b [ft 11^ |0> C2 |0> r2 + t?, |0> fl2 |1> C2 10>„ + Xl |0> fl2 |0> C2 11>,,] |0 fll Cl n > 



(21) 
(22) 



By tracing out the degrees of freedom of the qubits associated with the cavities ci , C2, and reservoirs f*i , r-i, the reduced density 
matrix of the bipartite two-qubit atomic system at time t (as determined by p, q) is obtained in the basis {|00) , |01) , |10) , |1 1)} as 



Pa u a 2 (t) = 



p ^ 

\b\ 2 (l -p) a*b(l-p) 

b*a (l-p) \a\ 2 (l-p) 

) 



(23) 



Likewise, by tracing out the degrees of freedom of the qubits associated with a\, «2, and r\, r2 (c\, Ci), the reduced density matrix 
of the remaining bipartite-equivalent partitions associated with the two-cavity (two-reservoir) system, p Cl , C2 (t) (P/-,,r 2 (0) can be 
obtained in the basis {|00> , |01> , |10) , |11)} by the substitution (1 -/>)-» q ((1 - p) -» yd)- 
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For the non-equivalent atom-cavity partition, we obtain an associated reduced density matrix (p a ,,c t ) (i=l,2) within the same 
subsystem as follows 



Pa u cSS) = 







o (n 



( \a\ 2 7d + b 2 

\ a\ 2 q \a\ 2 JT^^q 

| fl | 2 VT^pW N 2 (!-p) 

) 



Likewise the reduced density matrix (p 0j , rj ) associated with the atom-reservoir partition is obtained as 

( \a\ 2 q + b 2 



Pa,,r ; (0 = 



\a\ 2 y d N 2 VWV^ 

o N 2 Vwvr^ M 2 (1-p) o 





while the reduced density matrix (p Ci , n ) associated with the cavity -reservoir partition is obtained as 



( \a\ 2 {\ -p) + b 2 Q\ 

\a\ 2 y d H 2 V?V% 

o M 2 V?Vr^ M 2 <7 

o o o o ; 



(24) 



(25) 



(26) 



The reduced matrices (p a i,r 2 ^ j) associated with the atom-reservoir partition across different subsystems is obtained as 

(\a\ 2 (q + y d ) + \b\ 2 (l-y d ) 0^ 



\ b\ 2 y d a*b^y d {\-p) 

ab* Vr</d - P) \a\Hl-p) 
j 



(27) 



The remaining inter-system partitions, p auC ., Per,) (i + /)> are expected to display entanglement dynamics and non-locality 
characteristics not distinctly different from those considered above, and therefore will be omitted from further consideration. 



m. FIDELITY MEASURE FOR TWO-QUBIT PARTITIONS 

The fidelity measure < F\p(t),p(t')] < 1 quantifies the distance between the initial state p{f) and the evolved state p(/') at a 
later time t' and is defined as 12911 



F[p{t),p(f')]= Tr 
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(28) 



and satisfies the symmetry property, F[p(t),p(t')] = F[p(t'),p(t)]. Based on this definition, the fidelity measure F between the 
initial state at t=Q (for which p, q, yd=0) and which terminate as the state in Eqs.d23Tl. (l24T >. (l25l l, (l26l i or ( f27l ) at a later time t is 
obtained as 



l-- 6 e-^A 2 c 



\-e-'±(% + {tA c ) 2 + 4tA c ) 



FliPauaii^PauOiW)) = 1 ~ P> 
FliPcucM'PcucM = l-q, 

F 3 (p rur2 (0),p rur2 (t)) = \-y d , 

F 4 (p ai , C1 (0), p auCl (0) = ( ^a\\-p) + ^jb 4 + a 2 b 2 y)j 

yja 4 (l - p) + ylb 4 + a 2 b 2 q\ , 



(29) 



^a A e- ,4 r{tA c + A) 2 + 4 Jb 2 - -a 2 b 2 e-%(8 + (tA c ) 2 + 4tA c ) 



*( 1 + T) +y\b* + y 6 a 2 e-'±t 2 A 2 c b 2 



F 6 (p cun (0),p cun (t)) = b 2 +a 2 (\-p) 



b 2 + a 2 e- 
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where the corresponding fidelities near the vicinity of the exceptional point are provided within the square brackets. 

The fidelities F, (i=l to 6) between various subsystems for A c <4, are shown in Fig. [Tj At the exceptional point (setting 
ti - V =1), A c =4, and minimum passage time f=l, p—1. We note that the fidelity of the cavity-cavity partition experiences a 
minimum at t=l, which appears to be a characteristic feature of the cavity-cavity subsystem acting as a quantum channel. The 
local minimum in the fidelity measure may be linked to the merging of two eigenvalues at the exceptional point |3fJ and a 
corresponding loss of distinguishability. As expected, fidelities of all other subsystem decrease with time, and in particular the 
loss in fidelities of the atom-atom and reservoir-reservoir subsystems are almost similar. The decrease in fidelity of the reservoir- 
reservoir subsystem without any revival, is consistent with the role of a dissipative sink as there is a one-directional information 
flow between subspaces of P and Q, and the dissipative measure, jd increases with time. This feature is expected to be present 
in photosynthetic sinks as well. 




FIG. 1: Fidelities F, (i =1..6) between various subsystem partitions as a function of dimensionless time / and decay rate X c . F\ = atom-atom 
partition, F 2 = cavity-cavity partition, F 3 = reservoir-reservoir partition, F 4 = atom-cavity partition, F 5 = atom-reservoir partition, F 6 = 
cavity-reservoir partition. The unit system adopted here and in all other figures is based on fr = V =1, with time t obtained as inverse of Q.q (at 
A c =0) 




FIG. 2: Fidelities F, (i =1..6) of various subsystem partitions as a function of dimensionless time t at the vicinity of the exceptional point. 
Solid lines {F u F 4 ), DotDashed (F 2 , F 5 ), Dashed (F 3 , F 6 ) 



IV. NON-MARKOVIANITY BASED ON THE FIDELITY MEASURE 

For a system density matrix undergoing Markovian evolution under completely positive, trace preserving dynamical maps <f>, 
the fidelity measure F[p(t),p(t + r)] between the initial state p(t) and the evolved state p(t + r) at a later time t + r. satisfies ^2§i\ 



F[p(t),p(t + T)] = F[<X>(*)p(0),<X>(*)p(r)] 
^F[p(t),p(t + T)]>F[p(0),p(T)]. 



(30) 
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The violation of the inequality in Eq. ( l30b was proposed as a signature of non-Markovian dynamics by Rajagopal et. al. 
Accordingly negative values of the following fidelity difference function 



G(t,r) = 



F[p(t),p(t + T)]-F[p(Q),p(T)] 
F[p(0),p(T)] 



(31) 



are correlated with non-Markovianity during the evolution dynamics of the quantum state. Using the two-qubit matrices which 
appear in Eqs.d23Tl, (|24| |, 1241 . (l25l l and the expressions for p and q in Eqs.dTSb. (I161 l. the fidelity difference function G(t, r) can be 
evaluated as analytical expressions of t, t and A c . Due to the lengthy expressions for these functions, we only provide numerical 
results which show subtle changes of G(t, t) with A c for the various subsystems. 

Fig. |3 illustrates the fidelity difference G(t, t) as a function of t and t for the reduced density matrix corresponding to the 
atom-reservoir subsystem (Eq.(l25ll) for increasing A c . The sensitivity of the two-qubit dynamics on t translates in a striking way 
to the non-Markovianity measure, with the qubit dynamics becoming Markovian beyond a critical damping Ac, and as expected 
depends on the initial correlation amplitude factors, a, b. At A c =0, the oscillatory state of information exchanges between 
the atomic qubit and reservoir is captured by the symmetric pattern of non-Markovian dynamics alternating with Markovian 
dynamics. The non-Markovianity decreases with higher values of the cavity decay rate A c , consistent with decreased population 
transfer between the atom-cavity subsystem. Specifically, we note the persistence of non-Markovianity at small t , t values, even 
at higher values of the cavity decay rate A c « 5. 

As expected, there is insignificant non-Markovian signatures seen in a atom-cavity partition within the same subsystem. In 
contrast to the lack of non-Markovianity in the atom-cavity partition , we noted richness in non-Markovian dynamics for the 
atom-reservoir (p au r 2 ^ j) partition across different subsystems, and similar to that shown in Fig. [3] In this case, the two-qubit 
dynamics becomes fully Markovian at a smaller critical damping if a small initial correlation amplitude factor, a is used. 

The fidelity difference G{t, r) is also seen to undergo notable changes (with t and t) in the case of the reduced density 
matrix corresponding to the cavity-cavity subsystem ( Eq.(l23l) with the substitution (1 - p) — > q), as shown in Fig. [4] Here, 
the non-Markovianity appear to be enhanced in some regions as the cavity decay rate A c is increased. In the vicinity of the 
exceptional point, we note a transition in which regions of non-Markovianity merge, and Markovian processes are eliminated. 
The exceptional point appear to signal the enhancement of non-Markovianity in selected regions (of time t and A c ) for the 
cavity-cavity partition. 



a„=o 



,l„=0.1 




.l„=l 



A„=5 





FIG. 3: Fidelity difference G(t, t) as a function of t and r for the reduced density matrix corresponding to the atom-reservoir subsystem 
fEq.ll25t). and increasing A c . a=b=A=. Non-Markovian trends were not significant for the intra-system atom-cavity and cavity-reservoir 
partitions. Negative values indicating non-Markovianity are shaded green, while regions which are shaded purple imply non-violation of the 
inequality OOl l and hence Markovian evolution dynamics. Non-Markovianity decreases with higher values of the cavity decay rate A c . 



5 
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h 3 



\ c =a 




A„=0.5 



^=2 








FIG. 4: Fidelity difference G(t, r) as a function of t and t for the reduced density matrix corresponding to the cavity-cavity subsystem ( Eq.J23t 
with the substitution (1 - p) — » q). a=b=-^. Colour codes follows as specified in Fig. [3] Unlike Fig. [3] non-Markovianity continue to persist 
for higher values of the cavity decay rate A c in some regions. 



9 



A. Non-Markovianity based on other distance measures 

It may be worthwhile to examine other physical quantities, which are functions of the dynamical map, that describe the open 
system evolution of the physical states. In an earlier work fl411 , the use of relative entropy difference to qualitatively capture 
the departure from the completely positive Markovian semigroup property of evolution was verified to assume negative values 
under an open system non-completely positive dynamics. We recall that the relative entropy ll38ll of two density matrices p and y, 
defined by: S (p||y) = Tr[p(ln p - In y)] is positive and vanishes if and only if p = y. Under completely positive, trace preserving 
dynamical maps <I>, the relative entropy obeys monotonicity property which yields a relation which is analogous to Eq.(f30b 

S[p(t)\\p(t + T)] = 5[O(f)p(0)||<D(r)p(r)] 

< S[p(0)\\p(r)] (32) 

The relative entropy difference S (?, r) given by 

_ , S[P(0)\\P(T)1-S[p(t)\\p(t + T)] 

S(t,T) = (33) 

S[p(0)\\p(r)] 

is necessarily positive for all quantum states p(t) evolving under completely positive Markovian dynamics and violation of the 
inequality in Eq.(f32b i.e., S(t,T) < 0, denotes the occurrence of non-Markovian dynamical processes in the quantum system 
under study. 

Similar to the relative entropy difference, the trace-distance D[p\,pi\ - \\\p\ - pi\\ CD] where ||A|| = Tr[ YAtA] between 
pi, p2, can be used to define the trace-distance difference 

D(f . = D[p(0),p(r)-D\p(t),p(t + T)] 
D[p(0),p(T)] 

and to identify violation of the monotonically contractive characteristic feature associated with divisible Markovian mapping on 
the operator space. We note that both the relative entropy and trace-distance decline during a complete positive Markovian pro- 
cess, unlike the fidelity which increases under the same conditions. As pointed out earlier in Ref. 11411 . negative values of relative 
entropy difference d33b and the fidelity difference OTb only imply that the time evolution is not a completely positive Markovian 
process, positive values of these same quantities does not imply the occurrence of a Markovian evolution. Accordingly, the 
negative values of relative entropy and fidelity differences serve as sufficient but not necessary signatures of non-Markovianity 
(completely positive as well as non-completely positive). In order to highlight the latter point, we have repeated the calculations 
of the fidelity difference G(t, t) (Eq.OTTl) as a function of t, r done earlier, by using the trace-distance difference D{t, t) (Eq.d34b 
for the cavity-cavity two-qubit partition ( Eq.(l23]) with the substitution (1 - p) — » q), and illustrated in Fig. [5] 

Comparing the results of Fig.|4]and Fig. [5] we note that non-Markovianity is optimized in the vicinity of the exceptional point 
in both cases, however there are subtle differences. The trace-distance difference measure yields enhanced features in that it is 
able to detect the presence of non-Markovianity at small t which remains undetected by the fidelity difference measure. Similar 
trends have been noted ITill with the relative entropy difference, which shows the dependence of non-Markovianity dynamics 
on the metric entity. It appears that regions of contractive quantum evolution vary according to the measure used to detect 
violations of Markovian dynamics. In general one can conclude (at least on the basis of the results obtained in this study) that 
these differences are marginal and there is overall agreement in the critical regions of non-Markovianity. In the next Section, 
we make several important observations based on the non-Markovian patterns of Figs. [3] H] and and results of non-locality 
computed using the CHSH-Bell inequality. 



V. VIOLATION OF THE CHSH-BELL INEQUALITY 

The violation of the Bell inequality |S| < 2 is quantified by the CHSH- Bell inequality function S of a two-qubit density 
matrix where & = M(a, b) - M(a, b') + M{a' , b) + M(a', b'). M(a, b) is the correlated results (+1) arising from the measurement 
of two qubits in directions a and b. The CHSH-Bell inequality is violated when £ exceeds 2, and the correlations is considered 
inaccessible by any classical means of information transfer, while for values less than 2, the local hidden-variable theory can 
satisfy the inequality. Here, we investigate the role that the parameter A c play in the violation of a Bell inequality, with view 
to seeking a link with the non-Markovianity measure evaluated for various partitions in the earlier section. & is evaluated as a 
function of the correlations averages and for a given density matrix 
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Xc=0 .l„=0.5 Ac=l 1^=2 it=3.9 




123456 123456 123456 123456 123456 



t t t t t 

FIG. 5: Trace difference D(t, t) (Eg.J34t as a function of t and t for the reduced density matrix corresponding to the cavity-cavity subsystem 
( Eq.J23t with the substitution (1 - p) — » q). a=b=-j=. Colour codes follows as specified in Fig. [3] Unlike Fig. [4] non-Markovianity appear to 
be present at small r, t that remains undetected by the fidelity difference measure. In any case, there is overall agreement with Fig. [4] in that 
non-Markovianity is optimized in the vicinity of the exceptional point. 



using simple relations ll39ll 

£ = Max{S!,S 2 ) (36) 

Si = 2 \/4|M23| 2 + ("11 - "22 - M33) 2 

£2 = 2 V2|« 23 | 2 

Figure [6] a,b,c,d shows the regions (shaded purple) in which & exceeds 2, for the atom-atom, cavity-cavity, reservoir-reservoir 
and inter-system atom-cavity two-qubit partitions. Interestingly, the intra-system atom-cavity, atom-reservoir and atom-cavity 
partitions did not display any violation of the CHSH-Bell inequality for the range of time t and A c considered in Figure [6] 
Figure [6] c shows that the reservoir-reservoir subsystem exhibits a larger region of violation of the CHSH-Bell inequality when 
compared to all other subsystems. 

We can make several important observations by comparing the non-Markovian features of Figs. [3] IH [3 and the CHSH-Bell 
inequality violation trends in Figure [6] a,b,c,d. It is evident that non-locality present in one subsystem (e.g reservoir-reservoir 
subsystem) appears in conjunction with non-Markovian dynamics in an adjacent subsystem (e.g atom-reservoir or the cavity- 
cavity partition) with respect to changes in time t and decay rate A c . The cavity-cavity partition, for instance has notable non- 
Markovian features (Figs. [HE}, y et displays non-locality only in a narrow range of t and A c . There is also mismatch between 
non-locality and non-Markovianity in the inter-system atom-cavity two-qubit partition, with Bell nonlocal regions dominant at 
large A c > I at which the system dynamics is noted to be Markovian. While these trends seem applicable to two-qubit partitions, 
it remains to be seen whether higher dimensional qubit partitions follow similar trends. Our final observations relates to the 
match between Bell non-locality and non-classical correlations for select partitions which we illustrate next. 



(a) (b) (c) (d) 




1 2 3 4 1 2 3 4 0123456 1 2 3456 

t t t t 



FIG. 6: Left to right: S (a) for the atom-atom, (b) cavity-cavity, (c) reservoir-reservoir and (d) inter-system atom-cavity two-qubit partitions 
as a function of time t and A c . a=b=-^. (a-c) and for (d). The remaining intra-system atom-cavity, atom-reservoir and atom-cavity 

partitions did display any violation of the CHSH-Bell inequality (using Eq.d36b. Regions (shaded purple) for which S exceeds 2 are Bell 
nonlocal. The white lines divides the fine interplay between the two competitive terms which appear in the Bell function of Eq.j36ll 
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A. Classical and non-classical correlations 



In this section, we evaluate the correlation measure known as the quantum discord, D for the two-qubit state in Eq.d23l> and its 
cavity-cavity and reservoir-reservoir counterpart matrices. It is well known that the quantum discord is more robust as it captures 
nonlocal correlations not present in the entanglement measure I27ll40ll4lll . consequently it is non vanishing in states which has 
zero entanglement. The evaluation of D involves lengthy optimization procedures and analytical expressions are known exist 
only in a few limiting cases l42H44ll . 

The quantum mutual information of a composite state p of two subsystems A and B is given by l27l l40l l4ltl lip) = Sip a) + 
S (fis) - S (p) for a density operator in "Ha ® 'Hb- Pa (Pb) is the reduced density matrix associated with A (B) and S (p,) (i=A,B) 
denotes the well known von Neumann entropy of the density operator p,, where S(p) = -tr(plogp). The mutual information 
can be written as 



I(p) = S(pB)-S(p\p A ) 



(37) 



where S (p\pa) — S(p) - S (pa) is the quantum conditional entropy. 

After a series of measurements on A, the post measurement conditional state in B is given by ps\k = j^i^k ® WpCTt ® Ib) 
where the probability pk = \x[p{Ylk®^B)\ and {11*} denote the one-dimensional projector indexed by the outcome k. A conditional 
entropy of the of the conditional state of subsystem B, ps\k, is derived based on the cumulative effect of the mutually exclusive 
measurements on A as S (p| {11*}) = 2* fkS (pB\k)- The measurement induced mutual information is therefore given by 



i( P \{n k }) = s(pB)-S(p\{nk}) 



(38) 



Eqs.(|37]i and ( 1381 appear with slight differences due to incorporation of measurements in the definition for the latter. 

The classical correlation measure based on optimal measurements is evaluated using Ca(p) = suppj j -^(pHII*}), and is the 
maximum information about one subsystem that can be obtained via measurements performed on the adjacent subsystem. The 
quantum discord which measures the non-classical correlations is given by the difference in I(p) and C A {p) 

D A (p)=I(p)-C A (p) (39) 

Following the analytical approach in earlier works j42Ui4ll , expressions for the C and Q associated with the atom-atom 
bipartite system in Eq.(l23l is obtained as 



C(p aua2 ) = H(\a\ 2 (\-p)) 



1 



y jl-4\a\\l-p)p J 



D(p aua2 ) = H(\b\\l -p))-H( l-p) 



1+ ^1-41^(1 -p)p J 



(40) 



(41) 



where the function H(x) = -jclog 2 x - (1 - x) log 2 (l - x). Analogous expressions for C(p c ,,c 2 ), D (p C i,c 2 ) and C(p ri ,r z ), Dffir^) 
associated respectively with the cavity-cavity and reservoir-reservoir density matrices at time t, are obtained by respective sub- 
stitutions (1 — p) — ¥ q and (1 — p) — > jd in Eq. d40b and ( |4TT >. At t =0 and a = -4=, C(p fl] , a ,)=C(pa 1 ,a,) = l while the correlations 
of the cavity-cavity and reservoir-reservoir partitions are zero. Figures for the classical and quantum correlations as function of 
time t for the analytical solutions (p, q) in Eqs.(Q3]l, ( fT6] l are shown below. 

Figures Q and [8] show the expected shift in classical C(p) and quantum correlations D(p) away from the entangled atom- 
atom partition to the reservoir-reservoir partition. There exist larger regions (as function of time t and A c ) with zero quantum 
correlations than classical correlations. The appearance of D(p) is most significant in the reservoir-reservoir two-qubit partition 
as the system nears the exceptional point A c — > 4. Comparison of Figure|6]a,b,c and Figures|7]and [8] show that for the considered 
subsystem partitions (atom-atom, cavity-cavity, reservoir-reservoir), Bell nonlocal regions can be matched with appearance of 
non-trivial classical and non-classical correlations. This is most noticeable in the reservoir-reservoir subsystem, and once again 
we reiterate that these results appear valid for the two-qubit partitioned subsystems. 



B. Phase space approach of quantum dissipation 

The current study has highlighted some salient features linking entities based on the abstract Hilbert space , and within the 
confines of the two-qubit partitions examined in this work. However, the validity of the results of dynamics in various subsystems 
is expected to be dependent on the form of the quantum master equation utilized in Eq. (0. Possible artifacts which arise from 
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FIG. 7: Classical correlation C(p) for the (a) atom-atom, (b) cavity-cavity, (c) reservoir-reservoir two-qubit partitions as a function of time t 
and A c . a=b=-^=. 




FIG. 8: Quantum correlation Dip) for the (a) atom-atom, (b) cavity-cavity, (c) reservoir-reservoir two-qubit partitions as a function of time t 
and A c . a=b=-^p. There is transfer of the quantum discord from the atom- atom subsystem to the reservoir- reservoir partition. 

such Markovian forms have been examined by Kohen et al. JUt], who analysed time evolving density operators transformed 
using several approaches including Redfield theories, the master equations of Agarwal l46ll and the semigroup theory of Lindblad 
J*!], into Wigner phase space distribution functions. The Agarwal-Redfield (AR) equations of motion were seen to violate the 
semigroup form of Lindblad, showing that for certain initial states, the Redfield theory can violate simple positivity making it 
difficult to assign distinctly physical status for such states. 

The Agarwal bath ll46ll which is translational invariant and able to reach the thermal equilibrium limit, yields density matrix 
positivity for many initial conditions. However this bath system presents pathological effects as was noted in the density matrix 
negativity for a range of initial conditions at low temperatures 114711 in a model system of one primary oscillator. These obser- 
vations suggest possible association between the density matrix positivity and conditions which lead to the thermal equilibrium 
in system-bath correlations. In this regard, there is need to examine these links with incorporation of detailed attributes of the 
dissipative sink inherent in the the qubit-cavity-reservoir system of Eq. |[T), in future investigations. 

Finally, it is not entirely clear as to the role played by the bipartite partitioning of the six-qubit multipartite system (Eq. (l2ll ) 
in the observed matching or non-matching of regions with noticeable non-locality, non-Markovianity and quantum correlations 
in the subsystems examined in this study. Future investigations should be aimed at higher dimensional density matrices (e.g. 
tripartite systems) and among all possible combinations of partitions, to see if the trends observed in two-qubit systems are 
indeed universal. This can be a formidable task, which may require a rigorous approach lfl3l l48ll involving the mathematical 
formulations of the underlying abstract Hilbert space, to help form a basis for understanding the links between the quantum 
entities such as non-locality and non-Markovianity. 

VI. APPLICATION TO LIGHT-HARVESTING SYSTEMS 

The results of this study has immediate relevance to photosynthetic energy transfer in light-harvesting systems such as the 
Fenna-Matthews-Olson complex with moderate exciton-phonon coupling, and operation under physiological conditions. There 
is obvious similarity between the qubit-cavity-reservoir system in Eq. (Q]i and a photosynthetic model that constitutes the donor 
and acceptor protein pigment complexes with finite energy difference and a third party with a continuous frequency spectrum 
acting as the phonon-dissipative sink, and into which energy from the acceptor dissipates with time. The inclusion of compre- 
hensive details incorporating the solvent relaxation rate as well as kinetic behaviors of the electron transfer mechanisms and 
dynamic bath effects in Eqs. ( TTOb and ( fTTT i will extend the applicability of the model in Eq. I© to the more realistic environment 
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of photosy nthetic systems . 

Earlier studies lfl7j|49ll have shown that attributes such as the non- Markovian interactions and quantum correlations contribute 
via a delicate interplay of quantum mechanical attributes and the environmental noise in bringing about an optimal performance, 
of the light harvesting complex. Non-Markovian effects were seen to reduce the transport efficiency while increasing the life- 
time of entanglement lfl7ll by disrupting the optimal balance of quantum and incoherent dynamics required for efficient energy 
transfer. In another study involving light harvesting systems Jsoh . non-Markovian processes were noted to be most pronounced 
in the reorganization en ergy regime where the transfer efficiency is the highest, and linked with the preservation of coherence 
along critical pathways IU8I1 . 

Results based on the exciton entanglement dynamics of the Fenna-Matthews-Olson (FMO) pigment-protein complex l25ll . 
indicate increased oscillations of entanglement in the non-Markovian regime, with implications for a link between non- 
Markovianity and large coherence times. Interestingly, an earlier work on optimizing energy transfer efficiency by Silbey 
and coworkers l5ltl . showed that the interplay of coherent dynamics and environmental noise resulted in the optimal energy 
transfer efficiency at an intermediate level for several variables in light-harvesting systems. In particular, the reorganization 
energy and the bath relaxation rate played critical roles in displaying a non-monotonic-dependence on energy transfer. Spatial 
correlation effects were noted [51] to optimize the energy transfer only during strong dissipation, highlighting differences in the 
roles played by spatial correlation and temporal correlations. In this regard, there is need to further examine the subtle links 
between non-Markovian dynamics and the effective system-bath coupling range, taking into consideration the site energy based 
electronic coupling correlations f52ll . Future investigations along these lines can employ the model developed through Eq. ([8]), 
( [Tol l and dTTb as currently, the information on the exact mechanism by which non-Markovian processes act to bring about fast 
energy transport still remains unclear in light-harvesting systems. 

The inclusion of quantum effects involving non-locality and non-classical correlations will provide greater insight to the 
recently proposed optical cavity quantum electrodynamics setup ll53ll to investigate energy transfer mechanism in biomolecules. 

J I 

There is also scope to examine similar quantum effects in a viable Grover-like search process [54] via the exciton trapping 
mechanism at dissipative sink sites in photosynthetic systems. The incorporation of an approach incorporating non-Markovianity 
and the violation of the CHSH-Bell inequality may assist in developing new perspectives in interpreting experimental results 
obtained via two-dimensional electronic spectroscopy [550, 4-wave mixing measurements ll56ll and the hybrid optical detector 
1I57I1 which implements quantum detector tomography and characterizes simultaneous wave and photon-number sensitivities. 
To this end, future experimental work should incorporate the joint extraction of non-Markovian and Bell non-local features in 
biochemical systems, pending further refinement in optical techniques. 

VII. DISCUSSION AND CONCLUSION 

In summary, we have studied the links between non-locality, non-Markovianity and quantum correlations of two initially 
correlated atomic qubits, each located in a single-mode leaky cavity and interacting with its own bosonic reservoir. The appear- 
ance of non-Markovianity is associated with a negative value for the the fidelity or trace-distance differences, and act only as 
qualitative signatures of deviation from a complete positve Markovian behavior. Results of the differences in non-Markovianity 
based on the two measures: fidelity and trace-distance shows its dependence on the distance metric, and regions of contractive 
quantum evolution vary according to the defined measure used to detect violations of Markovian dynamics. 

The non-Markovian features in different two-qubit partitions (cavity-cavity and atom-reservoir partitions of the same subsys- 
tem, cavity-reservoir partition across different subsystem) show varying levels of erasure or enhancement of non-Markovianity 
depending on the cavity decay rate, A c . In particular, the exceptional point at which two eigenvalues merge, appear to signal the 
enhancement of non-Markovianity in selected regions (of time t and A c ) for the cavity-cavity partition. The changes in the vicin- 
ity of the exceptional point, also noted in quantum measurements l58ll . is significant as it may assist in distinguishing between 
processes with a classical origin and those with an intrinsically quantum mechanical coherence origin 15911 . 

Results of the CHSH-Bell inequality computed for various two-qubit partitions show that enhanced non-locality present 
in a specific subsystem (e.g reservoir-reservoir matrix) appears in conjunction with non-Markovian evolution in an adjacent 
subsystem (e.g atom-reservoir matrix). The mismatch between non-locality and non-Markovianity for a given partition is in 
contrast to the matching between non-locality and quantum correlations for regions spanned by time t and the cavity decay rate, 
A c . These results appear to provide weak signatures of subtle links between non-Markovianity, non-locality and non-classical 
correlations, bearing in mind that attempts to establish such links are very much reliant on the rigid formulation of the hidden 
variable model via the CHSH-Bell inequality relations in Eq. (f36b . Moreover the criteria used in this study provides only 
sufficient and not necessary tests of non-Markovianity, and also very much dependent on the initial state of correlations of the 
various subsystems. Currently, it is still not clear as to the origins of non-Markovianity as there is uncertainty as to whether 
information back-flow occurs due to the initial system-environment correlations or the intrinsic presence of specific dynamical 
features from which non-Markovianity emerge. Moreover, investigations involving tripartite states and higher dimensional states 
is needed to seek a definite conclusion on the perceived correlation/anti-correlation between non-locality, non-Markovianity and 
quantum discord. Lastly, the results obtained in this work highlight the need to adopt a unified model that incorporates non- 
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Markovianity, non-locality and non-classical correlations, as an effective approach in examining high efficiencies of energy 
transfer observed in light-harvesting systems. 
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